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PREFACE 


The purpose of this contract is to provide data processing and 
analysis in support of two 6SFC programs: tornado detection by radio 
frequency techniques and short pulse radar measurements of ocean waves. 
The first quarter's effort was concentrated primarily on the former 
program - specifically on the strip chart recording of the tornado 
detection raw data and preliminary analyses of the data. The next 
quarter's work will concentrate on digitizing, displaying, and analyzing 
selected data portions of interest based on the strip chart recordings. 



I. Introduction 


The tornado data analog tapes are 14 track tapes containing the 
information shown in table I. The recordings cover 4 HF and VHF 
frequency bands with vertical and horizontal polarizations and with 
linear and logarithmic scales. In addition the outputs from 2 
different lightning stroke detectors are recorded. Over 50 tapes 
have been currently generated by another contractor to GSFC. Recording 
speeds of 30 ips and 60 ips were used. 

Techno-Sciences has provided analog strip charts of all available 
tapes on GSFC supplied equipment in coordination with the contract 
technical officer. A complete log of this processing is included 
with this report in table II. Because the s’-. ip chart recorder has only 
7 useable channels, only half of the tracks on each tape were actually 
processed in most cases. The other tracks could have been processed 
by making a second pass through each tape. However, it was found that 
7 tracks had most of the useful information anyway. In particular, 
the logarithmic scale and timing tracks were found to be of limited 
usefulness. 

Current work is concentrating on digitizing selected analog tape 
portions for detailed analysis and display. 
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TABLE I 
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TORNADO ANALOG TAPE TRACK ASSIGNMENTS 


Track 

Mode 

Assignment 

1 

Di rect 

V-3 Mhz Lin 

2 

Di rect 

V-3 Mhz Log 

3 

Direct 

V-30 Mhz Lin 

4 

Direct 

V-30 Mhz Log 

5 

Direct 

V-VHF 139 Mhz Lin 

6 

Direct 

V-UHF 295 Mhz Lin 

7 

Direct 

H - VHF 139 Mhz Lin 

8 

Direct 

H - UHF 295 Mhz Lin 

9 

FM 

10 Khz Time Mark or 100 Khz at 60 ips 

10 

FM 

WWVB BCD Time Code 

11 

Di rect 

V-Lightning Stroke Det-Lin 

12 

Direct 

V-Lightning Stroke Det-Log 

13 

FM 

Taylor Tornado Det-Far 

14 

FM 

Taylor Tornado Det-Near 



II. Analog Strip C h art Kecordings 


3 


Table II is a processing log which sumriarizes the useful analog 
strip chart recording results by Techno-Sciences on GSFC supplied 
equipment. In addition to these results, other recordings were made 
but not kept because of various equipment problems which resulted in 
invalid or useless data displays. Many of the tapes were found to be 
very quiet, containing indiscernable signals or noise of a "white" 
rather than of the desired impulsive nature characterizing lightning 
discharges. Others, however contained considerable amounts of the 
desired impulsive noise across all frequencies. A short segment of 
one of these displays appears in figure "one". 

The data was recorded at a real time (one-sided) bandwidth of 
300 khz. As shown in Appendix B, the strip chart recording technique 
used is physically limited to be able to show only about 20% of the 
actual bandwidth. Even this could only be achieved at 1-7/8 ips 
playback speeds and 50 cm/sec strip chart speeds, which would generate 
unacceptably large amounts of output Hence, the strip chart recordings 
can only be viewed as a "quick look" tool. Selected portions of the 
data must be digitized for processing and display. 

As shown in Appendix C, the amount of digital data generated from 
even a short segment of analog tape is very large. Hence, the segments 
must be carefully selected. The addition of a timing track may be 
necessary in this regard as an aid in registration at some point. 
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III. Data Digitization 

One of the more interesting tapes (8/26/75 - #3) was selected for 
digitization and preliminary analysis and display. A request was 
made to GSFC for its digitization in September. However, as a result 
of the inavailability of equipment and appropriate software, the 
digitization has net been accomplished to date. Therefore an alternate 
digitization plan is being pursued. 

The Ampex PR-2200 tape recorder used in the Georgia Institute 
of Technology research instrun^ntation is temporarily available to 
playback selected tapes for digitization. This recorder is being 
used as an input to the A/D hardware developed at GSFC for data 
compression. Computer program^ are being written by Techno-Sciences 
to ingest the sampled data to core memory. Initial output will be 
to digital tape which has a basic (burst) limitation of 20,000 bytes 
per second (10 or 20 thousand samples per second dependent on quantiza- 
tion). Next programs will be developed to support rates to 100,000 
bytes per second for transfer to disk storage. 



IV. Plans for the Next Quarter 

Digitization will be accomplished as discussed in section III. 

As the digital data becomes available, analysis efforts will begin. 
The first analytical effort will be the isolation, display and 
analysis of individual lightning pulses. Time and frequency display 
will be accomplished using Fourier transforms. The time signal will 
be compared with earl i*?r published results. The transformed data 
will be investigated to determine actual bandwidth utilization. If 
the lightning stroke generates true impulses, the frequency content 
should be approximately flat to 300 khz and an individual pulse in 
tin« should approximate the impulse response of the system filter. 

Burst counts, inte-^ulse statistics and amplitude probability 
distributions will ba generated as an initial approach to determining 
a statistical model for lightning discharges. If sufficient data 
is available, tornado and non-tornado statistics will be segregated 
and contrasted for a possible tornado detection methodology. 
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Appendix A 

The Sampling Theorem 

Suppose that x(t) is an analog waveform containing frequencies 
from zero up to some maximum value B hz . The sampling theorem states 
that you must sample at a rate of at least 2B samples per second to 
reconstruct the waveform without error. The sampling 
rate, 2B, is called the Nyquist rate. To see the necessity for this, 
consider the simple co-sinusoids cos 2 (B-e )tv€nd cos 2 (B+e )t 
for any e . The nth sample at the sampling rate 2B per second of 
each is 

cos 2 IT (B- £ ) ^3 ■ (cos n ir ) cos + (sin n ir ) sin ^ ^ 

~ ( “1 / cos g 

and cos 2 it (B+ e ) = (cos n ir ) cos -(sin n it )s1n 

- (-1) cos — g — 

It is clear that there Is no way of distinguishing between the 
two waveforms based upon the sample values at the rate 2B. '^he 
characteristic problem of confusing two such waveforms due to lender- 
sampling is called al lasing or folding and can be avoided only by 
sampling fast enough or by filtering (with concurrent loss of high 
frequency content) to remove all frequencies higher than 1/2 the sampling 
rate. 

Conversely it can be intuitively justified that sampling at 
twice the highest frequency present is sufficient to reconstruct the 



Appendix A 
Page 2 

the waveform by arguing that the 2 samples per cycle can be used to 
find the 2 degrees of freedwn in a sinusoid - the amplitude and the 
phase. 

Quantitative proofs of the sampling theorem can be found in 
many textbooks on linear systems. 


Appendix B 


Strip Chart Recorder Display Limitations 

The Varlan strip chart recorder used for the processing of the 
analog tape recordings under this contract Is limited In usefulness by 
sampling theorem considerations presented In Appendix A. The plotting 
density of the strip chart recorder Is 200 points per inch, which by 
the sampling theorem allows for a maximum frequency resolution o' 
a waveform of 100 cycles per inch. The maximum strip chart speed Is 
50 cm/sec or approximately 20 Inches/sec. Thus the maximum frequency 
content of a signal which can be resolved Is approximately 2000 hz. 

At the recorded speed, the anlog tapes have frequencies up to 
300 khz • The maximum record speed used was 60 Ips. At a minimum 
playback speed of 1-7/8 ips, the highest frequency present Is then 
300/32 « 10 khz or 5 times the strip chart recorder's frequency 
resolving capability. 

Thus digital methods must be used to accurately represent and 
display the analog tape data. From the preceding considerations it 
Is seen that a minimum sampling rate of 20,000 samples per second 
must be used. Higher sampling rates are desirable to avoid analog 
tape flutter/wow problems at 1-7/8 ips by running higher playback 
(and hence higher sampling) speeds. 



Appendix C 


Digital Tape Limitations 

Based on the sampling rate considerations presented in 
Appendices A and B, it can be easily shown that great selectivity 
must be exercized in choosing the data to be digitized for display/ 
analysis. Each of the 14 300 khz analog tape tracks requires 
600,000 samples per second in real time or, at a recording speed 
of 60 ips, 10,000 samples per inch of analog tape per channel for a 
total of 140,000 samples per inch for all channels. Assuming 8 bit 
samples and continuous recoiding (no record gaps), a 2400 ft. 800 bpi 
digital tape can hold only (2400 x 12 x 800)/140,000 = 165 inches 
0 '^ analog tape or less than 3 seconds of real time sampled data. Hence 
the channels chosen for sampling and the actual segments of time 
should be carefully chosen for maximum information content. 



